Avalanche mixing of granular solids in a rotating 2D drum: 

diffusion and fractionality 

S. N. DorogovtsevQ 

A.F.Ioffe Physico- Technical Institute, 194021 St. Petersburg, Russia 
The dynamics of the avalanche mixing in a slowly rotated 2D upright 

drum is studied in the situation where the difference 5 between the angle of 
marginal stability and the angle of repose of the granular material is finite. An 
analytical solution of the problem is found for a half filled drum, that is the 
most interesting case. The mixing is described by a simple linear difference 
equation. We show that the mixing looks like linear diffusion of fractions un- 
der consideration with the diffusion coefficient vanishing when 5 is an integer 
part of 7T. The characteristic mixing time tends to infinity in these points. A 
full dependence of the mixing time on 5 is calculated and predictions for an 

experiment are made. 
PACS numbers: 64.75, 46.10 



The avalanche mixing problemGTQ, may be, is the only one among all numerous problems 
concerning the evolution of granular solids in rotating drumsiHi! that may be formulated 
with absolute mathematical clearness^!. 

Let us define first the so called avalanche mixing^. A partially filled 2D upright drum 
rotates slowly (see Fig. [I]). One can assume its radius be equal to 1, the angle of rotation of 
the drum plays the role of time. The free surface slope varies periodically between the angle 
of repose and the angle of marginal stability of the material. Small granules of different 
fractions (they differ only by color) move as a whole with a drum with one exception: 
periodically, when the free surface slope achieves the marginal stability value, the granules 
from the sector of the angle 5 near the upper half of a free surface fall down (see Fig. [T] - 
granules fall down from sector A to sector B). Here 5 is the difference between the angle 
of marginal stability and the angle of repose. Therefore, the only reason for mixing in the 
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problem is granular flow (avalanches) between the wedges. One should describe, how the 
mixing develops in time. 

We assume that the fractions becomes be mixed homogeneously after the fall - the 



appears to be near the experimentally realized situational. When 5 — > 0, the mixing can be 
described in the frames of a simple geometrical approach0~i, and an analytical theoryi'S can 
be used. There is no such clearness, when 5 is finite!'!. Below we consider analytically the 
finite S case in the most interesting situation of a half filled drum. (When a drum is half 
filled and 5 = 0, the mixing is absent - see papersEHi, so in this situation, the effect of 
nonzero 5 is most striking.) 



FIG. 1. The avalanche mixing scheme in the case of a half filled drum. When the free surface 
slope attains the marginal stability angle value, the granules of different fractions flow quickly 
from sector A to sector B, undergoing mixing. Then the free surface slope turns to be equal to the 
repose angle of the material under consideration. The angle 5 is the difference between the angle 
of marginal stability and the angle of repose. The mixing fractions are not shown. 

Let us assume that there are two fractions - black and white. Then the mixing dynamics 
is described by a set {c m }. Here c m is the concentration of the black fraction in sector B 
(Fig. [I]) after the m-th avalanche falls, that will be at the moment t = mn/S. In the case 
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of a half filled drum, the equation for c m appears to be very simple: 

C m = jl - (j ~ | ) } C m ^ [n/5 ] +{■$- | } c m-[7r/<5]-l (1) 

(it can be understood easily from Fig. |l[). Here [•] denotes an integer part of a number. The 
values c m , m = — [ir/5], — [ir/5] + 1, . . . , — 1, are used as an initial condition. We shall use 
the notation n = [tt/S] for the integer part and a = n/5 — [n/5] for the fractional part. 

It is not a simple task to write directly an explicit form of the solution of Eq. (|lj), if 
n > 2. Let us find the solution in the case where 5 <C it, i.e. n ^> 1. The idea of the trick 
to be used is obvious from the following form of Eq. ([!]) : 

Cm+n C m d{c m C m _i) . (2) 

One is tempted to introduce different variables for the different scales: a la time variable t 
for a large scale in the left-hand side of Eq. (0) and an angle variable p for a short scale in 
the right-hand side. 

We shall pass on carefully to continual variables (the necessary notions will be made 
later). First, let m be continual now. Then ip = mS. We shall define the variable t in 
the following way. Let us calculate the following difference: c m+n — c m _ g ( a ), where q(a) is 
unknown function for the moment. 

S(q(a) - a)^ + j {(1 - a)g 2 (a) + a(l - g(a)) 2 } . (3) 

Now we can choose q(a) demanding the coefficient of the first derivative be zero. Therefore, 
q(a) = a and one may write 

dc dc 

c m+n -c m _ a ^(n5 + a5)— T = 7c— T . (4) 

at at 

Hear the time variable t is introduced in the following way. We observe periodically the 
system under consideration with the periodicity interval 7r (like using the stroboscopic effect) 
and describe the result of the observation by the continual variable t. 



Finally, we get the partial differential equation 

dc(<p,i) 5 2 a(l — a) d 2 c(<p,t) 



(5) 



dt 2ir dcp 2 

One can easily verify that the boundary conditions for Eq. (||) are c((p = 0) = c(cp = tt) and 
dc((p = 0)/d(p = dc(ip = n)/d(p. The initial distribution of the fractions defines the initial 
condition for Eq. ([|). Thus, the avalanche mixing is described by linear diffusion equation 
(H) with the diffusion coefficient that vanishes at the angles 5 = ir/k, where k is integer. 
Note, that previously, the resonance-like behavior of some other quantity due to the same 
reason was found numerically!'!. 

The solution c((p, t) of Eq. (|jD solves the problem, since, of course, one can find trivially 
all c m from c(ip,t). We shall not present this expression here but shall write immediately 
the answer for the mixing starting from the following initial configuration. Let it be a thin 
layer of the black fraction over the white fraction in the drum at the initial moment. Then, 
in our situation, one may use the initial condition c((p,t = 0) = ^5{ip). Here 5((p) is the 
5-function and 47 defines the amount of the black fraction. Then the solution is of the form 



z{<p, i) = 4 7 I - + - exp \ -k 2 - 



5 2 a(l — a)t \ cos (2k(p) 



[l eJ^-t 2 a(l-a)t) , (6) 



where 9 3 is the theta-function. If t <C 7r/ (25 2 a(l — a)), 

5sj2a(l-a)t I 25 2 a{\ - a)t/u ) 

For longer time one obtains exponential relaxation with the relaxation time r (i.e. the 
mixing time): 



r- 1 = 2 -5 2 ( n - 

71 \0 
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(8) 



to a homogeneous mixture Cqo = 47/71". (Of course, expression (|3|) is valid for any initial 
conditions.) 
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Note that continualization of difference equations is very dangerous trick (e.g. see the 
discussion of the topic in the handbook of E. Pinneylll). Nevertheless, in our specific case, it 
works! We shall not try to justify it rigorously, (one of the reasons of its availability is fast 
decreasing of all derivatives of c((p, t) and the other one is the separation of the scales - see 
Eq. (fj)) but shall compare the main derived result, i.e. the mixing time, with that obtained 
in a more direct way. 

According to usual prescriptions (e.g. see handbook^) , one may search for the solution 
of Eq. ([!]) or Eq. (||D in the form of a linear combination of the terms A™, where Xj are the 
roots of the characteristic polynomial 

A n+1 - (1 - a)X - a = . (9) 




FIG. 2. Scaled inverse avalanche mixing time (tt/25 2 )t~ 1 vs. tt/5. The solid line is obtained 
from the "exact" solution of the characteristic polynomial The dashed line represents expres- 
sion (||) which is, formally speaking, obtained in the limit of small 5 <C vr. 

If is the multiplicity of the Xj root, then the corresponding contributions to the general 
solution of Eq. (|2|) are m S: > \Xj\ m cos(u jm), m S] |Aj| m sin(u;,,m), where Sj — 0, 1, . . . , fij — 1 
and Uj = arg(Xj). Thus, to obtain the long time relaxation we should find the maximal 
|Aj| = Xmax which is lower then the root A = 1 of Eq. (§). The mixing time is expressed in 
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the terms of \ ma x'- t^ 1 = (1/5) \og(l / X max ) , so it may be calculated directly for any value of 
a and a finite n. The result is shown by the solid line in Fig. 0. The dashed line represents 
the answer (|8|) obtained for small angles 5 <C n. The closeness of these curves (even for large 
angles S) looks impressing. It verifies the quality of the method used above. 

In an experiment, the angle 5 usually fluctuates with time. One can see easily that if 
the deviations from the mean value of 5 are greater then a very small angle 6 2 /n {S/tt <C 1), 
then they should smear the dependence t(5) and eliminate the singularities. It is not an 
easy task to take into account exactly these fluctuations. The most naive estimation for 
the averaged value looks like 5 2 /3ti. In the experiment^, 5 = 8° ± 2°, so from our 

estimation, r ~ 80 turns of the drum. Unfortunately, we did not find any experimental data 
on this quantity for a half filled drum, so the comparison is impossible at the moment. 

In conclusion, we have described analytically the dynamics of the avalanche mixing in 
the case of a half filled drum where the difference between the angle of the marginal stability 
and the angle of repose plays the principal role. The mixing looks like diffusion with the 
diffusion coefficient vanishing when 5 is an integer part of tt. 
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